Correlated atomic motion, phonons and spin-ordering in bcc ^He 

N. Gov and E. Polturak 

Physics Department, 

Technion- Israel Institute of Technology, 

Haifa 32000, Israel 






o 

in 



> 

m 

o 

-(— • 



o 



% 



We propose a new way to treat nuclear magnetism of solid He. We argue that the magnetic 
interaction arises indirectly as a consequence of correlated zero-point motion of the ions. This 
motion lowers the energy of the ground state, and results in a coherent state of oscillating electric 
dipoles. Distortion of the electronic wavefunctions leads to hyperfine magnetic interactions with the 
nuclear spin. Our model describes both the phonon spectra and the nuclear magnetic ordering of bcc 
^He using a single parameter, the dipolar interaction energy Eq. The model yields correctly both 
the u2d2 symmetry of the ordered phase and the volume dependence of the magnetic interaction. 



PACS: 67.80.-s,67.80.Jd,67.80.Cx 

The spin-ordered phase of bcc ^He presents a 
ficult challenge to accurate theoretical descriptiont 
The main uroblem is to explain why the transition 
temperatureB of 10~^K is two orders of magnitude larger 
than the nuclear dipolar interaction ~10~^K. The preva- 
lent description is in terms of atomic exchange cy- 
cles that produce competing ferro and antiferromagnetic 
interactionsu. This model has conceptual problems that 
arise from the fact that in order to fit experimental 
data, many exchange cycles involving a large number of 
atoms are needecL|Jn fact, it is not sure that this ex- 
pansion convergeajfl. Additionally, an overall consistent 
descripti|Qn of the experimental data has still not been 
achievedtliJ. It is therefore of interest to consider another 
approach. We propose that magnetic ordering is due to 
correlations in the zero point atomic motion, which how- 
ever does not involve exchange of atomair. These correla- 
tions lower the energy of the ground stateH and result in a 
coherent state of oscillating electric dipoles, which modify 
the phonon spectrum. Finally, this motion produces an 
oscillating magnetic polarization of the electronic cloud 
which interacts with the nuclear spin. We show below 
that this hyper-fine type interaction has the right order 
of magnitude, and leads naturally to the distinct u2d2 
antiferromagnetic phase. This scenario is somewhat sim- 
ilar to the mechanism of hyperfine enhaBced nuclear mag- 
netism that occurs in certain insulatorsB. The calculated 
volume dependence of the magnetic interaction is in good 
agreement with experimental data. 

At temperatures which are high compared to the mag- 
netic interactions (T^lmK), the local motion in bcc '^He 
can be treated in the same way as in '*Hd3. The usual 
Helium inter-atomic potential is the Van-der Waals in- 
teraction, due to the zero-point fluctuations of the s- 
shell electrons. This interaction is calculated assuming 
that the nuclei are stationary (Born-Oppenheimer ap- 
proximation), which is reasonable since the nuclear mo- 
tion is usually negligible. The result of averaging over 
random electric dipolar fluctuations leads to an isotropic 
1/r^ attraction. In addition there is the hard-core repul- 



sion due to the overlap of the filled electronic shells. This 
interaction, strictly valid for an isotropic medium, is also 
used in anisotropic cases, such as for atoms in a crystal. 
We would like to describe the lowest correction to the 
isotropic interaction resulting from the zero point mo- 
tion of the He ions. This motion, resulting from confine- 
ment, is especially large and anisotropic in solid He. In 
the bcc phase, this motion is large along the major axes 
(100,010^pi), due to a shallow, double-minimum lattice 
potentialEnl3. If we relax the Born-Oppenheimer approx- 
imation, and allow some relative motion between the ion 
and the electrons, we will obtain that the zero-point mo- 
tion of the ion creates an oscillating electric dipole. One 
possible source of this relative motion may be the large 
difference in inertia between the ion and electrons. If 
this ionic motion is random and isotropic it would have 
no observable consequences. However, we will show that 
the solid phase can gain energy if this ionic motion is 
correlated. 

Directional oscillation of the nucleus will break the ro- 
tational symmetry of the nuclear position relative to the 
electronic cloud. This is equivalent to a mixing xif the 
s and p electronic levels. The size of this mixingLl is of 
the order of the square root of the ratio of the oscil- 
lation frequency of the nucleus compared with the fre- 
quency of the s — > p electronic transition, and comes out 
~ 10^"". The distortion of the electronic wavefunction 
is described through the mixing coefficient A, defined by 
|f/;) ~ |s) -|- A |p).We therefore have in the bcc phase an 
isotropic Van-der Waals interaction with an additional 
dipolar interaction acting along the major axes. If the 
relative phases of the oscillatory motion between differ- 
ent ions are correlated the net interaction energy can be 
non-zero, and becomes negative for special 'antiferroelec- 
tric' configurations, in which the phase of the oscillations 
of neighboring atoms differs by tt (see Fig.l). This con- 
figuration has a zero total dipole moment, as required. 
The dipolar potential for the case, of correlated motion 
has a long range, 1/r^ dependence^. The ground-state of 
the crystal in which the zero-point motion of the atoms is 



correlated may be describgd as a global state of quantum 
resonance (coherent statelll) between the two degenerate 
configurations shown in Fig.l, eacj^ of which minimizes 
the dipolar interaction energy EiUp 
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where Tq is the position of an atom and the summation 
Yi runs over the entire lattice. The atomic electric dipole 
moment |^| induced by the coherent motion is propor- 
tional to A through 
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We would like to describe the effect of this addi- 
tional interaction on the phonon spectrum. Except 
for the Ti(llO) branch, the phonons of the solid are 
well describedJsy the Self-Consistent Harmonic (SCH) 
approximatioriij. The additional part is given by an ef- 
fective HamiltonianLl which describes the coherent dipo- 
lar interactions. This Hamiltonian was used to de- 
scribe the, interaction between excitons in a dielectric 
mediumOtj, where the flipping of a dipole out of t 
ground state is treated as a local mode with bare energ; 
Eq = 2 \Edip\ ■ The excitations of this Hamiltonian are a 
periodic modulation of the relative phases of the dipoles, 
with respect to the ground state array (Fig.l). Due to 
the lower symmetry of the dipolar array compared to the 
crystal symmetry a coupling between the usual phonons 
and the |dipolar modulation exists only along the (110) 
directionlZI. The hybridization of the dipolar modulation 
and the usual phonons results in a soft Ti(llO) mode 
with approximately half the energy given by SCH calcu- 
lation. At the edge of the Brilloin zone this phonon has 
the energy Eq . In addition there appears a localized exci- 
tation of energy 2Eq, which is invoWed in mass diffusion 
and contributes to the specific-heatQ. The excitations in 
our picture are therefore three phonon branches (only 
one of which is renormalized) and a localized excitation 
which is a quantum analogue of a point defect. 

In the case of bcc '''He the bare local mode energy Eq 
was measured directly by NMR of a ^ He impurity, and 
the phonon spectrum was measured using neutron scat- 
teriiiff. Our calculation is consistent with both sets of 
datcu. There are no neutron scattering measurements of 
the phonon spectrum of bcc '^He, but there are sound 
velocity data that indicate that the slope of thfi T i(llO) 
phonon is about half of the SCH calculatioEllil^ll3. Since 
we predict that this ratio shoulijndeed be 0.5, we take 
half the energy of the calculateclij SCH Ti(llO) phonon 
at the edge of the Brillouin zone to be the bare dipole- 
flip energy Eq ~ 5K at V=21. Bern's/mole. According 
to our model, the energy of the localized mode involved 
in thermally activated self-diffusion is 2Eq ~ lOK. This 
value is in excellent agreement with the activation en- 
ergy measur|ed by x-ray diffraction, ultrasonics and NMR 
experimental^ at V=21.5cm^/mole. A similar activation 



energy is also obtained from, the excess specific heatll3, 
and pressure measurementala. The localized excitation 
of energy 2Eq is further described in Ref . [7] . 

We now describe the magnetic interaction arising from 
our model. The lowest \p) level of the He atom has the 
electrons iiua spin S ^ 1 state due to strong exchange 
interactiontil, of the order of 0.25eV. In addition this level 
is spht into 3 sublevels with J — L^ S — 0, 1, 2. The 
splitting, due to spin-orbit coupling, is of the order of 
~ 1.5K. In the ground state, the \p) electrons will reside 
in the ^P2 sublevel, with an oscillating magnetic moment 
Me of size ^ XrUe, where m-e is the magnetic moment 
of an electron. Because part of the magnetic moment 
is now in the \p) state, there appears a net uncancelled 
moment of equal size in the \s) component of the elec- 
tronic wavefunction. In ^He the nuclear spin is zero, and 
the magnetic moment of the electrons has no effect. In 
contrast, in "^He the nuclear magnetic spin / — 1/2 will 
interact with the oscillating electronic magnetic moment, 
mainly due to the contact term of the \s) electron at the 
nucleus. The magnetic interaction is of the hyper-fine 
typelli, and the energy associated with it, Emag, is given 
by 



E„ 



Stt 



Me • MnSir) 



(3) 



where M„ is the nuclear magnetic moment and the 
calculation of the matrix element is done just as for the 
\s) level of the hydrogen atomtll. We show below that 
the maximum value of Emag/kB — 0.75mK (for V=24 
em's /mole), jj^^,]^ larger than the direct nuclear dipole- 
dipole interaction, and is of the right magnitude to ex- 
plain the high transition temperature of nuclear ordering 
in bcc "SHc. 

We would like to make a quantitative calculation of 
the strength of the magnetic interaction as a function of 
pressure at T=0. In our model, the magnetic energy (|^) 
will change with pressure due to changes of the electronic 
magnetic polarization Mg. First, the pressure changes 
the admixture of the \p) orbital in the ground state wave- 
function \ip), through changes of the mixing parameter 
A. Second, as the solid is compressed, the 3 sub-levels 
(J — 0,1,2) broaden into partially overlapping bands. 
This overlap results in a reduction of the net alignment 
of the electronic magnetic moment Mg. 

The mixing parameter A depends on pressure through 
y/VEo, where V is the molar volume. This behaviour 
arises from the definition of Edip (|]),(0) and the factor 
|ro — r^l in this equation is proportional to the molar 
volume V . We ..take Edip as a function of pressure from 
the measuredllja activation energy 2Eq. A is found to 
increase with decreasing molar volume V . The effect of 
the broadening of the sub-levels with pressure can be 
approximated using the overlap integral of the three sub- 
levels 
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where the overlap integral F{V) = J ipl[E)'il) -i[E)dE . 
Here, i/jHE) and -0-1 (^) are taken as normalized gaus- 
sians centered at the energies of the spin-orbit psJigned 
and anti-aligned levels, which are ^ 1.5K apartEJ. The 
effective energy width 7 of ipi{E) and ip-i{E) is esti- 
mated using simple band calculation, as the change of 
the Coulomb energy of the \p) electrons due to-pverlap- 
ing \p) wavefunctions on neighboring He nucleillj 
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where ipo ('') j i'i (^) are the p-state wavefunctions of the 
electrons around the central atom and its i neighbors. 
In this calculation we also take into account the large 
spatial spread in the atomic position inside the wide po- 
tential wellEI, of approximately ±0.8A. We find that 7 
ranges from ^OK for V=24 cm^/mole to 5.5K for V=19 
cm'^/mole. The overlap factor 1 — F{V) is very sensi- 
tive to volume, changing from 1 at V=24 cm^/mole to 
~ 0.01 at V=19 cm^/mole. One can see that as the 
volume decreases the broadening of the bands increases 
thereby decreasing the net magnetic polarization of the 
oscillating electronic cloud. 

The strength of the magnetic interaction should be 
proportional to the measured Curie- Weiss temperature 
9. In Fig. 2 wc compare the magnetic spJitting 2E„iag (0) 
with the two available sets of measuredEl 6. We find that 
the volume dependence of Emag agrees very well with 
that of 6, i.e. 2Emag = const ■ 9. The constant is 1 for 
the data in Fig. 2a, and 2 for the data in Fig. 2b. 

We now consider the symmetry of the ordered spin 
system. The existence of the hyper-fine splitting means 
that the simple quantum resonance condition on each 
site is broken. Now, the two antiferoelectric configura- 
tions shown in Fig. 1 are not degenerate, with an energy 
difference of 2Emag per site. It is possible to restore the 
degeneracy of the overall ground state, and hence the 
quantum resonance condition. This can be done if the 
nuclear spins become ordered as well, in a spatial config- 
uration which has the same interaction energy with each 
of the two states of Fig. 1. The electric dipolar inter- 
actions in the ground-state of Fig.l are such that each 
simple cubic sublattice of the bcc has no net interaction 
with the other sublattice. Two possible static arrange- 
ments of the nuclear spins that fulfill the resonance re- 
quirement on each sublattice are shown in Fig. 3. Both 
arrangements are the simplest which ensure that there 
is an equal number of atoms with electronic and nuclear 
spins aligned (and anti-aligned) in both degenerate con- 
figurations of the electric dipoles. These arrangements 
preserve the overall time-reversal symmetry of the sys- 
tem at zero field. Since these configurations exist on both 
sublattices, we end up with an u2d2 arrangemepts which 
is the symmetry of the ordered nuclear phaseE3IHJ. Spon- 
taneous symmetry breaking due to ordering of the spins 
will select one of the major axes (100,010,001), similar 
to a ferro or antiferromagnet. The electric dipoles oscil- 



lating along the two orthogonal axes will have no mag- 
netic interaction energy, and are therefore in the usual 
quantum resonance of Fig.l. We point out that in our 
model, the u2d2 symmetry of the spin ordered nuclear 
phase is a natural consequence of the coherent-state of 
the oscillating dipoles (Fig.l), and does not depend on 
the magnitude of the magnetic interaction. 

To conclude, our model enables us to describe both 
the phonon spectra and the nuclear magnetic ordering of 
bcc ^He using a single parameter, which is the thermal 
activation energy Eq. Experimental values of E^ were 
measured in several experiments. The model describes 
correctly both the symmetry of the ordered phase and 
the volume dependence of the magnetic interaction. 
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FIG. 1. The two degenerate 'antiferroelectric' dipole ar- 
rangement in the ground-state of the bcc phase. The arrows 
show the instantaneous direction of the dipoles. 




FIG. 3. The two simplest u2d2 static nuclear spin (arrows) 
arrangements that maintain the quantum resonance of both 
the electric and magnetic dipoles (Fig.l). Only the spins on 
one simple-cubic sublattice of the bcc are shown, with the 
other sublattice having identical arrangement. 
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FIG. 2. The magnetic splitting 2Emag (Eq.3) as a func- 
tion of volume, calculated using different experimental data 
sets for the activation energy Eo: dashed line [18], dash-dot 
[17] and solid line [16]. The solid circles in (a) and (b) are 
the experimental Curie- Weiss temperature 6 [1]. To compare 
with the higher set of data (b) , we multiplied our calculated 
splitting by an extra factor of 2. 



